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Summary 

Multiple genetic loci have been associated with body mass index (BMI) and obesity. The aim of this study was to 
investigate the effects of established adult BMI and childhood obesity loci in a Greek adolescent cohort. For this purpose, 
34 variants were selected for investigation in 707 (55.9% females) adolescents of Greek origin aged 13.42 ± 0.88 years. 
Cumulative effects of variants were assessed by calculating a genetic risk score (GRS-34) for each subject. Variants at the 
FTO, TMEM18, FAIM2, RBJ, ZNF608 and QPCTL loci yielded nominal evidence for association with BMI and/or 
overweight risk (p < 0.05). Variants at TFAP2B and NEGRI loci showed nominal association (p < 0.05) with BMI 
and/ or overweight risk in males and females respectively. Even though we did not detect any genome- wide significant 
associations, 27 out of 34 variants yielded directionally consistent effects with those reported by large-scale meta-analyses 
(binomial sign p = 0.0008). The GRS-34 was associated with both BMI (beta = 0.17 kg/m 2 /allele; p < 0.001) and 
overweight risk (OR = 1.09/allele; 95% CI: 1.04— 1.16; p = 0.001). In conclusion, we replicate associations of established 
BMI and childhood obesity variants in a Greek adolescent cohort and confirm directionally consistent effects for most of 
them. 
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Introduction 

Prevalence of overweight and obesity in adolescents has in- 
creased over the last decades, and the observed rates in 
Greece are higher or comparable to those reported for other 
European countries (Lobstein & Frelut, 2003; Tzotzas et al., 
2008). Although the increased obesity levels have been at- 
tributed to environmental changes, a strong genetic com- 
ponent has also been shown to contribute. Genome-wide 
association studies (GWAS) have been successful in identify- 
ing multiple genetic loci associated with BMI and/or obesity 
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(Frayling et al, 2007; Scuten et al., 2007; Loos et al, 2008; 
Thorleifsson et al., 2009; Wilier et al., 2009; Speliotes et al, 
2010; Bradfield et al, 2012). 

Large GWAS meta-analyses, from the Genetic Investiga- 
tion of Anthropometric Traits (GIANT) (Speliotes et al., 
2010) and Early Growth Genetics (EGG) consortia (Bradfield 
et al., 2012), have identified multiple BMI-associated loci. 
The GIANT meta-analysis (Speliotes et al, 2010) of 249,796 
individuals confirmed 14 established obesity susceptibility 
loci and identified 18 new loci associated with adult BMI. 
Twenty-three of the 32 variants also showed direction- 
ally consistent effects on children's BMI. The EGG meta- 
analysis replicated associations for seven established adult 
BMI variants and identified two novel loci robustly associ- 
ated with increased childhood obesity risk (Bradfield et al., 
2012). 
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In this study, we examined the effects of established 
adult BMI and childhood obesity associated loci reported 
by GIANT (Speliotes et al, 2010) and EGG (Bradfield 
et al., 2012) meta-analyses respectively on BMI and over- 
weight risk in Greek adolescents. In total 34 single 
nucleotide polymorphisms (SNPs) were selected for inves- 
tigation; 32 were based on the GIANT loci (FTO, MEM18, 
MC4R, GNPDA2, BDNF, NEGRI, SH2B1, ETV5, 
MTCH2, KCTD15, SEC16B, TFAP2B, FAIM2, NRXN3, 
RBJ, GPRC5BC, MAP2K5, QPCTL, TNNI3K, SLC39A8, 
FLJ35779, LRRN6C, TMEM160, FANCE, CADM2, 
PRKD1, ERP1B, PTBP2, MTIF3, ZNF608, RPE21A, 
NUDT3) (Speliotes et al., 2010); and two were based on 
the EGG loci (OEFM4, HOXB5) (Bradfield et al, 2012). 



Materials and Methods 

Our sample comprised 707 (55.9% females) adolescents of 
Greek origin aged 13.42 ± 0.88 years. Participants were 
drawn from the TEENAGE (TEENs of Attica: Genes and 
Environment) study. A random sample of 857 adolescent stu- 
dents (Table SI) attending public secondary schools located 
in the wider Athens area of Attica in Greece were recruited 
in the study from 2008 to 2010. Details of recruitment and 
data collection have been described elsewhere (Ntalla et al., in 
press). Prior to recruitment all study participants gave their ver- 
bal assent along with their parents'/guardians' written consent 
forms. The study was approved by the Institutional Review 
Board of Harokopio University and the Greek Ministry of 
Education, Lifelong Learning and Religious Affairs. 

DNA samples were genotyped using Illumina HumanOm- 
niExpress BeadChips (Illumina, San Diego, CA, USA) 
at the Wellcome Trust Sanger Institute, Hinxton, UK. 
Genotype calling algorithm used was Illuminus (Teo et 
al., 2007). Sample exclusion criteria included: (i) sam- 
ple call rate < 95%; (ii) samples with sex discrepan- 
cies; (iii) samples with genome-wide heterozygosity of 
<32% or >35%; (iv) duplicated/related samples identified 
by calculating the genome-wide pair-wise identity by de- 
scent (IBD) for each sample using PLINK (Purcell et al., 
2007); from each pair with a tt a > 0.2 the sample with 
the lower call rate was excluded; (v) samples with evi- 
dence of non-European descent as assessed by perform- 
ing multidimensional scaling (MDS) in PLINK (Purcell 
et al, 2007) by combining the TEENAGE dataset with 1184 
individuals from the HapMap phase III populations. SNP ex- 
clusion criteria included: (i) Hardy Weinberg Equilibrium 
(HWE) exact p < 0.0001, (ii) MAF < 1%, (iii) call rate < 
95% (for SNPs with MAF > 5%) or call rate < 99% for 
SNPs with MAF < 5%. Genotypes were imputed using the 
directly typed data and phased HapMap II genotype data from 
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the 60 CEU HapMap founders using the program IMPUTE 
(Marchmi & Howie, 2010). 

Body weight (kg) was measured to the nearest 0.1 kg with 
the subjects barefoot and dressed in light clothing by the use 
of a weighing scale (Seca Alpha, Hamburg, Germany). Height 
was measured to the nearest 0. 1 cm using a portable stadiome- 
ter with participants being barefoot with their shoulders in a 
relaxed position, their arms hanging freely and their head in 
a normal position, with the eyes looking straight ahead. BMI 
was calculated as weight (kg)/height squared (m ). Subjects 
were classified as normal weight, overweight and obese ac- 
cording to the age- and sex-specific criteria adopted by the 
International Obesity Task Force (IOTF) (Cole et al., 2000). 

Before testing for associations, BMI was log-transformed to 
achieve normal distribution. The association of each variant 
with BMI was tested using linear regression. Overweight and 
obese subjects were classified into one category, and over- 
weight risk was tested using logistic regression. Obesity risk 
was also tested with the use of logistic regression. All analyses 
were adjusted for age and sex and were carried out assum- 
ing an additive genetic model in SNPTEST (Marchini et 
al., 2007). In order to investigate whether some of the se- 
lected variants have sex-specific effects, stratification analyses 
by sex were also performed. In order to investigate the cu- 
mulative effects of variants, a genetic risk score, comprising 
all selected variants for investigation (GRS-34) was calculated 
for each subject by adding the effect alleles of the SNPs. 
For imputed variants, we used best-guess genotypes to calcu- 
late the score, using a genotype probability threshold of 0.8. 
Given that it has been proposed that allele weighting may 
have a limited effect (Janssens et al., 2007), we did not weight 
risk alleles for their individual effect sizes. Linear and logis- 
tic regression models were performed in STATA-11 (Stata- 
Corp, College Station, TX) to investigate the effects of the 
score on BMI and overweight risk respectively. Quanto vl.2.4 
(http://hydra.usc.edu/gxe/) was used for power calculations. 
A two-sided p < 0.05 was used as the threshold for nominal 
significance, andp < 5 x 10~ 8 was used as the threshold for 
genome-wide significance. 



Results 

Three variants corresponding to previously identified GIANT 
loci {FTO, TMEM18, and FAIM2) yielded nominal evidence 
for association with both BMI (beta ± SE: 0.57 ± 0.19, p 
= 0.001; beta ± SE: 0.46 ± 0.24, p = 0.031; and beta ± 
SE: 0.49 ± 0.20, p — 0.015 respectively) and overweight risk 
(OR = 1.33, 95% CI: 1.06-1.67, p = 0.019; OR = 1.46, 
95% CI: 1.08-1.97, p = 0.011; and OR = 1.29, 95% CI: 
1.03-1.63, p = 0.025 respectively) in TEENAGE (Table 1). 
Variation at the QPCTL and ZNF608 loci yielded nominal 
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evidence for association with BMI (beta ± SE: 0.54 ± 0.25, 
p = 0.026 and beta ± SE: 0.38 ± 0.25, p = 0.047 respec- 
tively), while the index variant at the RBJ locus was associated 
with overweight risk (OR = 1.32, 95% CI: 1.05-1.66, p = 
0.017). The strongest association was observed at the FTO 
locus with BMI, which also accounted for the largest propor- 
tion of variation (0.96%) (Table 1). Obesity risk was nominally 
associated with variation at the FAIM2 locus (OR = 1.58, 
95% CI: 1.05-2.40, p = 0.025), and with a variant at the 
BDNF locus (OR = 1.72, 95% CI: 1.00-2.95, p = 0.028) 
(Table S2). Despite the lack of statistically significant evidence 
for association for the majority of variants, overall 27 of the 34 
GIANT and EGG loci yielded directionally consistent effects 
in TEENAGE (binomial sign test p = 0.0008) (Table 1). 

In the sex-stratified analysis, BMI and obesity were also 
nominally associated with a variant at the TFAP2B locus (beta 
± SE: 1.22 ± 0.39, p = 0.002; OR = 1.96, 95% CI: 1.08- 
3.55, p — 0.026) (Table S3a) in males. In females, a variant 
at the NEGRI locus also yielded nominal evidence of asso- 
ciation with risk of obesity (OR = 3.15, 95% CI: 1.10-9.02, 
p = 0.018). Variants at MTCH2, LRRN6C and TMEM160 
were also nominally associated with overweight or obesity 
risk. However, the direction of effects observed in our study 
for these variants was inconsistent with those reported from 
the GIANT consortium (Speliotes et al., 2010) (Table S3b). 

The GRS-34, which was constructed to investigate the 
cumulative effects of the 34 variants under study, was sig- 
nificantly associated with BMI (beta = 0.17 kg/m 2 /allele; 
p < 0.001) and explained 3.2% of BMI variation. The differ- 
ence in mean BMI between subjects with the higher GRS-34 
(> 37 effect alleles) (2.0% of subjects) and those with the lower 
GRS-34 (< 22 effect alleles) (2.3% of subjects) was 4.3 kg/m 2 
(Fig. 1). Consistent with the observation for BMI, the asso- 
ciation of the GRS-34 and risk of overweight showed that 
each additional risk allele was associated with an 1.09-fold in- 
creased odds of overweight (95% CI: 1.04-1.16; p = 0.001). 



Discussion 

Our findings regarding the association of established adult 
BMI associated loci with BMI and/or risk of overweight 
in adolescents reflect previous reports either on single pedi- 
atric cohorts (den Hoed et al., 2010; Bradfield et al., 2012; 
Zhao et al., 2011) or large consortial meta-analyses (Bradfield 
et al., 2012). In the EYHS study, most of the BMI associated 
loci identified by GWAS in adults were also associated with 
anthropometric traits in children (den Hoed et al., 2010). In 
a case-control study of European American children, Zhao 
et al. (Zhao et al., 2011) reported evidence for association 
with common childhood obesity for nine of the 32 GIANT 
loci, while 28 of the 32 loci yielded consistent directional 
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Figure 1 Distribution of the genetic risk score and cumulative 
effects of the risk alleles from 32 GIANT BMI and 2 EGG 
childhood obesity SNPs. Mean (±SE) values for BMI (kg/m 2 ) 
are also presented. 

effects. The EGG consortium meta-analysis (Bradfield et al., 
2012) also replicated robust associations for seven established 
adult BMI loci (FTO, TMEM18, POMC, MC4R, FAIM2, 
TNNI3K and SEC16B) with risk of childhood obesity; and 
identified two novel childhood obesity loci both of which had 
yielded directionally consistent effects in the GIANT meta- 
analysis of adult BMI (Speliotes et al., 2010). Consistency of 
directionality for the great majority of tested variants indicates 
that a large genetic component of BMI and obesity overlaps 
in children and adults (Bradfield et al., 2012). 

Our study has a small sample size compared to both con- 
sortial meta-analyses, and power to confirm associations even 
at nominal significance level for all tested loci is relatively low. 
The effect sizes at the GIANT-identified loci require sample 
sizes that range from a few thousands to a few hundred thou- 
sand in order to achieve >80% power to detect association. 

Despite the small effect sizes of established adult BMI asso- 
ciated variants, it has been shown that they have cumulative 
effects on BMI and obesity risk (Li et al., 2010). In our study, 
the combined effects of the 34 adult BMI and childhood 
obesity variants reported by the GIANT and EGG meta- 
analyses respectively were significantly associated with both 
BMI and risk of overweight, and the effect sizes are compara- 
ble to those reported in large meta-analyses of adults (Speliotes 
et al., 2010). 

In this report, we investigated the individual and cumula- 
tive effects of established adult BMI and childhood obesity 
associated loci with BMI and overweight risk in a sample of 
Greek adolescents. We report evidence of nominal association 
for several loci and show that cumulatively tested variants are 
associated with both BMI and overweight risk. Our findings 
also support evidence for a large shared genetic component 
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between adult and childhood BMI and obesity and validate 
the TEENAGE study as a cohort in which to study the ge- 
netics of anthropometric traits. 
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